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AB ST RAC T 


The number of zeros contained in a given interval (a,b) for a 
solution to a Sturm-Liouville differential equation is of importance in 
many problems of mathematical physics. This number maybe deter¬ 
mined through Sturm’s Comparison Theorem, Given one zero of the 
solution to a Sturm-Liouville differential equation, a technique, based 
upon Sturm’s Theorem, of computing the next consecutive zero of the 
solution is proposed. The existence of a function which satisfies the 
desired end results of the proposed technique is shown. The technique 
is then applied to Bessel’s differential equation and the results tabula¬ 
ted for the first 20 roots of and Jj(x), Unfortunately this 

technique did not achieve the desired result of convergence to succes¬ 
sive zeros of the given Bessel Function. 

The writer wishes to express his appreciation for the assistance 
and encouragement given him by Professor E. J. Stewart of the U. S, 
Naval Postgraduate School in this investigation. 
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1. Introduction. 




Consider a differential equation system. 

[r(x) y' (x)] + P(x) y(x) = 0 (1. 1) 

y(a) = 0 y(b) = 0 

where r(x) and p{x) are positive, continuous, real valued functions 
on [a,b]. The system (1. 1) is known as a Sturm-LiouvilfeSystem. 

A frequent problem encountered in mathematical physics requires 

;V'' 

the determination of consecutive characteristic roots of a system (1. 1)l 
The solutions to (1. 1) usually do not possess the nice regular charac¬ 
teristics of the trigonometric functions; consequently the former are 
in general more difficult to manipulate. 

When the value of one characteristic root a of (1. 1) is known, 
Sturm's Comparison Theorem suggests a possible means of computing 
the next consecutive root b in the following sense: 

A sequence of functions is generated from (1. 1). Each element of 
the sequence is obtained by selecting a particular value of the indepen¬ 
dent variable x = and substituting this value into the coefficient P(x) 
and r(x) of (1.1) giving 

[r(xpy'(x)]' + P(x.) y (x) = 0 (1.2), 

The solution of (1. 2) is forced to vanish at y = a, and the next consec¬ 
utive zero of this solution is chosen as next element of the 

sequence. It is desired that the sequence so generated result in an 
equation of the form (1. 2) such that its solution vanishes at x = a 
and X = b, the zeros of the solution of (1. 1). 

The theory is developed for a general Sturm-LiouvilfeSystem but 
specific applications are limited to Bessel's Equation of order n . 

(n a positive integer). 
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We shall define an indexing of the roots of a solution to a differen¬ 


tial equation as follows: 

Let rj be the first root > 0, and r^ the next root greater than 
r j. Furthermore, if r^ is the i th root ^ 0, then is the 

(i+l)st root > 0. 

If is the given root of a solution to a differential equation, we 

shall define r, , , as the next consecutive root of the solution, 
k+1 


2 


< 8 # 4i » * I -MH* <1 ^ I mm —H >1 

«ili«i ,1 |t MM( 

<iP *1*1 iitl**^ |i> 

tadill 

nn •! 4? « 



I 



2. Mathematical Background, 

This section contains theorems and mathematical developments 

which will be called upon in later sections. 

2, 1. Sturm's Comparison Theorem. 

Given two equations of a Sturm-Liouville type: 

y'' + Ay=0 (2.1.1) 

y(a) = 0 

u" + Bu = 0 (2. 1. 2) 

u(a) = 0 

where A and B are positive, real constants and B is greater than A. 

Let us determine a solution y(x) to (2. 1. 1) such that y(a) = 0. 

The elementary theory of differential equations tells us that the general 
solution is of the form 

y(x) = C j sin ^Ax cos "^A x . 

When X = a 

y(a) = Cj sin cos-yTa = 0 

C, sin'^L a 

C 2 = ■pzT' 

cos *yA a 

Since Cj^ is an arbitrary constant, choose cos VTa. 

y(x) = k^ [sin '\fAx. cos Va a - C OS VTx sin VTa] 
y(x) = sin VX (x - a). 

Similarly, (2. 1. 2) will have a solution u(x) = kj sin Vb(x - a) . 

The assumption A < B tells us u(x) will "oscillate" more rapidly than 
y(x) or that u(x) will have at least as many zeros as y(x) on a given 
interval (a, b). 

When the coefficients of (2. 1. 1) are positive, real valued functions 

of the independent variable x, the solutions obtained must be defined 
in terms of so-called "special" functions such as Bessel functions. 
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Charles Sturm (1803-1855) made a study of the rates of "oscilla¬ 
tion*' of solutions to equation#of the form 
[ r (x) y'(x) ]' + P (x) y(x) = 0 
[ r (x) w'(x) ]' + Q(x) w{x) = 0 
y(a) = w(a) = y(b) = 0 

where r(x)^0; P(x) and Q(x) are continuous on [a,b]. He proved 
that, in general, the larger P(x) the faster y(x) will oscillate. More 
concisely, Sturms Comparison Theorem is stated as follows: 

Theorem 2. 1 Sturm's Comparison Theorem. 

Given two equations; 

[ r(x) y']' + p(x) y = 0 (2.1.3) 

and [ r(x) v' ]' + q(x) v = 0 (2. 1.4) 

where r(x) is positive in the closed interval [a,b]; 
r(x), p(x) and q(x) are continuous in [a,b] and 

q(x) > p(x) with strict inequality for at least one point in 

[a,b]. If a solution y(x) of (2. 1. 3) has consecutive 

zeros at x = x and x = x, , with x, > x , and x, , x 
o 1 1 o 1 o 

are in the interval [a,b], then a solution v(x) of (2. 1.4), 
which vanishes at x = x^, will vanish again in the open 
interval (x^, Xj) . 

Proof: 

Suppose y(x) > 0 in (x^, Xj) and y'(x^)>0, y'(Xj) < 0 and 
v'(x^)>0. Then if y(x) and v(x) are solutions to (2. 1.3) and 
(2. 1.4) respectively, we have the identities 

[ r(x) y'(x) ]' + p(x) y(x) = 0 (2. 1. 5) 

and [ r(x) v'(x)]' + q(x) v(x) = 0 (2. 1.6) 
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If we multiply equation (2. 1. 5) by -v(x) and equation (2. 1.6) by 
y(x) we have 

-[r(x) y'(x) ]' V (x) — p(x) y(x) v(x) = 0 (2, 1. 7) 

[r(x) v'(x)]' y(x) + q(x) v(x) y(x) = 0 (2. 1.8) 

Adding (2. 1. 7) and (2. 1. 8) then integrating from x^ to x^ we 
have 


|Xi Xj 


•(x)[y<x)v'(x) - y'(x)v(x)] | +/" [q(x)-p(x)]y<x)v(x)dx = 0 

o J-x. 


or 


rx, 

r(Xj)y'(xj)v(xj) s / [q(x)-p(x)] y(x) v(x) dx (2. 1. 


9) 


o 

since y(x^) = 0 if x^ is a solution to (2. 1.3). 

Suppose now that v( x ) > 0 everywhere in (x^, Xj), (i. e. if 
v(x^) = 0 and v* (x^) > 0 then v(x) has not vanished anywhere in 
(x^, Xj) ), then this implies that r(xj) v(xj) y'(xj) is negative while 


X 

(q(x) • p(x) ) y(x) v(x) dx is positive. This contradiction implies 
o 

that v(x) cannot be positive everywhere in (x^, x^). Similarlyi the 
assumption v'(x^) < 0 and v(x) negative everywhere on (a,b) again 
leads to a contradiction and the theorem is established. 

2.2 The Normal Form. 

Consider a differential equation 

y" + P(x) y' + Q(x) y* 0 (2. 2. 1). 

Let y = u V, then y* s uv' +u'v 



and 


y‘% uv'* + 2u'v' + u"v. 

Substituting lor y in (2.2.1) we have 

uv" + [2u' + P(x)u]v' + [u" + P(x) u' + Q(x)u] v s 0 (2.2.2). 
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Suppose we choose u such that the coefficient of v* vanishes. Set 


2u* + P(x) u = 0 
du/u + dx = 0 

u = C exp [- 

Let C = 1 and we have 
u* = ^ P{x) u 

u” = ^P{x) u* - ^ u ^^^dx 

^ (P(x))^ 1 dP/j 

= + — - / dx u . 

Substituting u' and u" in (2.2. 2) we get 


uv 


" + t-P(x) u + P(x) u] v' + [Q(x)u + u - i ^^/dx u 


P(x)' 


u ] V = 0 . 


We see that we have a factor of u 
divided out leaving 


in every term,therefore it may be 


v" + [Q{x) - ^ -^ ] V = 0 (2. 2. 3) 

We say that (2. 2. 3) is equation (2. 2. 1) transformed into normal 
form. We will see that the normal form plays an important role in our 
later development. 

2. 3. Existence Theorem. 

Consider a differential equation reduced to normal form 


y" + R(t)y=0. (2.3.1) 

y(a) = y{b) = 0. 

We desire to know if there exists a t in the interval (a,b) such 
that the equation 

u" + R(F) y = 0 (2. 3.2) 

R(F) is R(t) evaluated at t = F such that a solution u{t) to (2. 3,2) 
can be forced to have zeros at both t = a and t = b. The answer to 
this question may be stated in the form of a theorem as follows: 
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Theorem 2. 2. 


Given a differential equation of the form: 

y** + R(t) y = 0 (2. 3.3) 

where R(t) is positive, real valued, continuous and monotone in (a,c), 
where b < c, and whose solution y(t) has consecutive zeros at t = a 
and t = b, a < b < c. Then there exists a t in the interval (a,b) such 
that the equation 

u*» + R(F) u = 0 (2. 3.4) 

has a solution u(t) with a zero at t = a and moreover, u(t) has its 
next consecutive zero at t = b. 

Proof: 

First assume R(t) is monotone decreasing, then 

R(a) ^ R(t) > R(b) for all t in (a,b). Choose that solution u^(t) of 

u** + R(b) u = 0 such that u^(t) vanishes at t = a. A solution that 

satisfies this requirement is ^(t) = sin y^R(b) (t-a). Clearly 

u^(a) = 0 and moreover, u^{t) will have its next consecutive zero, 

for t greater than a, when the argument of the sine function is equal tt, 

that is, when VFb) (t~a) = ir . Since R is assumed decreasing, 

R(b) < R(t) for all t in (a,b[therefore, Sturm*s Comparison Theorem 

tells us that will vanish again at t = t^ where b < t j. Similarly, 

let u** + R(a) u = 0 have a solution u (x) which vanishes for t = a. 

a 

Then with R(t) < R(a) for all t in the interval (a,b), w^ have that 
Ua(t) will vanish again at some point t = t^ such that a < t^ < b. 

Case II. 

Now suppose R(t) is increasing. Choose u^(t) a solution of 
V** + R(b) V = 0 

such that a root at t = a. R(t) < R(b) for all t in (a,b), hence v^(t) 
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will vanish again when t = < b. Similarly, v" + R(a) v = 0 has a 

solution V (t) which vanishes at t = a and will have its next consecutive 

ci 

zero at some t = t., b < t., 

4 4 

In either case, R(t) increasing or R(t) decreasing, the root b of 
equation (2. 3. 3) is between two roots of (2. 3. 4) determined by evalua¬ 
ting R(t) at t = a and at t = b. This and the continuity of R(t) assures 
us that there exists a t in (a,b) such that ut^ (t), a solution to (2. 3.4) 
which vanishes at t = a, will vanish again when t = b, and the theorem 
is established. 

An interesting conclusion from Theorem 2. 2 is the following: 
Theorem 2.2. 1. 

Given a system of the type (2. 3. 1) where R(t) is monotone, then 
the next consecutive root of a solution u{t) to (2, 3.4) such that 
u(a) = 0 must lie in the interval (a + — - - , a = ) . 

y R(b) yRla) 

The proof follows directly from the proof of Theorem 2. 2. The end 

points of the interval (a + ^/yR(b), a+ '/r(I) ) are precisely the 

t. found in Theorem 2. 2. 

1 

Note that the arbitrary constants occuring in the solutions to equations 
(2. 3. 1) and (2. 3.4) have no effect upon the zeros and, therefore, will 
be ignored in all further considerations. 

We should also note that in a general Sturm-Liouvilfe System (2. 3. 1), 
the function R(t) is not always monotone. 
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y”+xy* + (x^- n^) y = 0 


3. Bessel’s Equation. 

Let us now consider a typical Sturm -Liouville System: 

(3. 1) 

y(a) =: y(b) = 0 (3. la) 

Equation (3. 1) is the general form of Bessel’s Differential Equation of 

order n. A solution of (3. 1) is 

y(x) = J^(x) (3. 2) 

where J (x) is known as the Bessel function of the first kind, 
n 

We see that before we can derive any benefit from the theorems of 
section 2, we must reduce (3. 1) to its normal form. Write (3. 1) in 
the form 

1 .2 

(3. 3) 


y'» + i y» + (1 _ IL ) y = 0. 

X ^ 

X 


Section 2. 2 tells us that (3. 3) may be reduced to normal form by 
the transformation 

y = V exp [ - 2 dx ] 

= V exp [ - 2 logg X ] 

(3,4) 


= V X 


-1 /2 


Then 


, 1 , -3/2 . , -1/2 

y = - T v(x) + v' X 


(3.5) 


y" = j V X 


-5/2 1 , -3/2 1 , -3/2 ,,-1/Z 

: - v’x - — v'x + v”x 


(3.6) 


Substituting (3.4), (3.5), and (3. 6) into (3.3) and simplifying we 


have the desired normal form of Bessel's Equation 

II + n + V _ Q 

V" + ( 1 + - y- ) V = 0. 

4x 


(3. 7) 


The transformation (3.4) implies that (3,7) has a solution 

v(x) = x^ J„(x). (3. 8) 


n 


Trivially any zero of (3, 2) is also a zero of (3„ 8). 
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Consider the differential equation 

+ w = 0 (3. 9) 

w(a) = 0 

We know that a solution to (3, 9) 
w(x) = k sin(x-a) 

vanishes at x = a and again at x = a + ir . 

Suppose we examine (3.7) and (3.9) with respect to Theorem 2. 1 
(Sturm's Comparison Theorem). Let us denote the coefficient of v in 
(3. 7) as 

P(x) = 1 + {3. 10) 

4x 

We see from (3. 10) that we must consider two cases, 1) when n = 0 
and 2) when n 1 . 

Case I. 

n = 0. We see that P(x) > 1 and continuous, therefore, theorem 
2. 1 applies. We conclude that if the solution (3. 8) to (3. 7) vanishes 
when X = a, it will vanish again on the interval (a, a + tt). This implies 
that the difference between consecutive roots of is less than ir . 

Moreover, we see that 
lim P(x) = 1 . 

X 0=3 

This tells us that with respect to the index i of roots r^ ; 
lim 1 I = IT . (3. 11) 

Case II. 

1 '>J Z 

n 1. In this case 0 < P(x) < 1 when ^ ^ri -1 <x< . 

Then x must be in the interval ( ^ V4n^ -1 , o® ) in order that (3, 7) 
have an oscillating solution. Theorem 2. 1 tells us that a solution to 
(3. 7) which vanishes at x = a will not vanish again in (a, a + ir ). 
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Therefore, the difference between the consecutive roots of J (x) is 

n' ' 

greater than tt . Similarly, 

lim P(x) = 1, when n ^ 1 , 

tells us that the limiting difference is the same as in Case I, namely the 
difference given in (3. 11)« 

The treatment of the differences between consecutive roots of the 
Bessel functions in the above two cases, suggests that further treatment 
of the problem also be made in two parts, i« e« 1) n = 0 and 2) n ^ 1 
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4. Application to Bessel's Equation of Order Zero« 

When n = 0, the normal form of Bessel's differential equation 
reduces to 

v" + (1 + — ^ ) V = 0, (4. 1) 

4x 

Suppose the conditions 

v(a) = v{b) = 0 (4. 2) 

are imposed. Let the coefficient of v in (4. 1) be denoted by P{x), i. e, 

P(x) = 1 + (4. 3) 

4x 

Clearly P(x) p> 1 for all x, P(x) is continuous and decreasing^ Then, 

from theorem 2. 2, there exists an x in the interval (a,b) such that 

the equation 

u" + (1 + —^ ) u = 0 

4x 

has a solution 

u (x) k sin + ~irz 

with a root at x == a and, moreover, will have the next consecutive 
root at X = b. 

Let us now attempt to construct a sequence which will converge 

to X « There are two immediate choices for the first element of the 

r I' 

sequence, a or a + tt « We will generate the sequence by means 

of the following recursion formulae: 


P(x ) - 1 + 

(4.4) 

4x. 


1 


u!* f P(x.) u. = 0 

1 'll 

(4,5) 

u.{x) = k sin \/^P(x^)(x-a) 

(4.6) 

x., , = a + it/Vp(x.) 

1+1 1 

(4, 7) 
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Wherein (4, 6) is a solution to (4. 5) which is forced to have a root at 
X = a; and is the next consecutive root of (4.6). We see from 

(4.4) that P(x^) ^ 1 for all x. . This implies that < a + ir for 

all x^ . 

Case I. 


Choose x^ = a, giving 


p(Xo) - 1 + —= 1 + ^ 

4x 4a 

o 

(4.8) 

u‘* + P(x ) u =0 
o o o 

(4.9) 

Uo(x)*= sinYP(x^)(x-a) 

(4. 10) 

Xj = a + ir /Y"P(x^) 

(4. 11) 

Since P(a) P(x) for all x in the interval 

(a,b), theorem (2.1) 

(Sturm's Comparison Theorem) tells us that 

Xj lies in the interval 

a < x^ < b . 


Then Xj is substituted in the recursion equations giving 

p(Xi) = 1 t ^ 

4xi 

(4. 12) 

u'l + P(xj^) Uj = 0 

(4. 13) 

Uj^(x) - kj sin yP(xj^) (x-a) 

(4. 14) 

x^ = a + IT / yP(x^ j 

(4. 15) 


We cannot apply theorem 2. 1 with respect to (4. 13) and (4.3) 
since we cannot guarantee the necessary inequality between P(Xj) and 
P(x). However, we may apply theorem 2.1 with respect to (4. 13) 
and (4. 9) since P(xj) < P(x^). This tells us that x^ ^ Xj . 
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If we expand \4„ 15) 


IT 


X2 = a. + 


% 


1 + 4x, 


= a+ j — 

Y ^ 


/ 

4/ a + 


^ 4a' 


/ 


we see the sequence is generating a rather strange appearing continued 
fraction; strange in the sense that the more common continued fractions 
do not contain radical expressions^ Moreover, note that the fractional 
form becomes exceedingly more complicated as each new term is added„ 
Let us now assume that for some m >0, we have 

*m-l > ’‘m-2 

The inequality (4„ 16) implies that the following inequality also holds 

(4. 17) 


P(x J = 1 + -V- < 2 ) 1 + 

m-1 . ^ m-^ 4 Y ^ 

m-1 m-^ 


Theorem 2« 1 immediately tells us that x is greater than x . 

^ m ® m-1 


Inductively we have shown that the sequence is increasing. 
Moreover, we have shown that the sequence is bounded above by a + tt, 
consequently, the sequence converges to some 3c such that 

TT 


X == a + 


1 + 


1 

4x 


t ' Tt ? »2 ^ c '.'’2 2 

(x-a) - 4x TT 


4x > 1 


{4x 4 l)(x - a)^ 


^ 2 2 y-w 

4x TT = 0 


4x ^ - 8a x ^ + {4a^ + 1 - 4 it^) x ^-2a 4. + a^i= 0 , 

Then x is the root of (4. 24) that is nearest a + -it 

^ 2, 


(4„23) 

{4o 24) 

We see from 


(4, 23) that for large roots, 4 x 1 and the large roots of (4o 24) may 
be approximated by: 
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4 ^*" ^ . .Z z z _ 

X (x-a)-4x ir=0 

/•- x2 2 n 

(x - a) - TT =0 

<ri^ O 

X = a + TT 

A simple calculation will show for any given root *a* of J^(x) that 
(4. Z4) has a root near a + tt . 

Case IL 

Let us choose x = a 4 tt , 
o 


P(x ) = 1 + 
o 


1 


4(a+ irr 


u*‘ + P(x ) u =0 
o o o 


u (x) = k sin '/pix T (x-a) 


(4.25) 

(4.26) 


Xj = a + IT / 




(4.2 7) 

(4. 2 8) 

We have P(a + it) < P(x) for all x in (a,b)j therefore, we know x 
cannot be in (a,b), therefore, we see 


1 


b < Xj < a + IT 

Substituting x^ in the recursion formulae we have 

1 


P(Xj) = 1 + 


4x 


1 


Uj" + P(Xj) u^ 


Uj(x) = kj sin '/p(x^) (x-a) 
x^ = a + -rr / V*^F*(x^) 


(4. 29) 

(4. 30) 

(4. 31) 
(4.32) 
(4. 33) 


The inequality (4.29) tells us that Pfxj) ^ P(x^), and vre get x^ < x^ 
from Theorem 2. 1. Moreover, (4. 29) also tells us that Pfxj) < P(x) 

for all X in (a,b) and x^ J> b . 

If we reverse the assumed inequality in (4. 16) we can show induc¬ 
tively that this sequence £xj;' , generated by starting at x^ - a + ir , 
is a monotone decreasing sequence. Moreover , we see that |x^| > b 
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for all i, hence is bounded below. Then the sequence converges to an 
X|^ where 

= a + TT / '/p(x^ ) 


tr 



(4 + l){xj-a)^ - 4 - 0 

4 - 8a Xj^ + (4a^ + 1 - 4tr^) - 2a - 0 (4« 34) « 

The quartic (4o 34) lias coefficients which .are identical to the coef¬ 
ficients of like powers of in (4,24) ^ therefore we conclude that 
Xj = X , This shows that there is no difference in the limit between 
the sequence generated in Case I and that generated in Case II, Un- 
fortunately, we have also shown that x > b , but we have no indica.tjOn 
of the difference between x and b. 

The quartic (4, 24) has been solved with the aid of the CDC 1604 
Digital Computer for the first 20 roots of J^(x), The computed values 
of X are compared with the tabulated values of the roots of « The 

difference between x and the tabulated value of b is shown in the error 

column. The difference has its maximum at the first root and decreases 
quite rapidly as x grows larger. 

Simple calculations show that ^ is almost midway between b and 

a + ir for the first root, but this relation does not persist as the root 
number increases^ In fact x approaches b more rapidly than a + tt , 

The fact that the sequence x, does not converge to b implies that 

some change to the recursion formulae must be made^ The most direct 

method of changing the formulae consists of a modification to P(>;)„ 
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4. L Linear Average* 


Suppose a sequence 



were to converge such that 


(4* 1. 1) holds, 


i. e. 


b 


a + TT 



where 


is the limit point of a sequence 
{A% ^ + l){b-a)^ = A% 

4% ^[(b-a)^ - ir^] = - (b-a) 



4 


(4. 1. 1) 



We have shown that (b-a) < tt , therefore, we know that % is real* 

Computed values of % are displayed in table II, along with the tabulated 

values of b. The error column is the difference between x and M , 

and indicates the amount x must be decreased if the sequence is to 
converge to b* Noting that % is almost midway between a and x the 

following modification to the recursion formulae is suggested: 



1 


P(x.) 

= 1 + 2 

(4. 1.2) 

1 

4(a + X.) 

■ 


- 1 - 


u. f 

1 

P(x.) u. = 0 

1 1 

(4. 1.3) 

u.(x) 

= sinV^(x.) (x-a) 

(4. 1.4) 

^i+1 

s a TT / 

(4. 1.5) 


In essence, the modification consists of merely substituting 
{a 4 X ) ^ 

' i' in place of x. in {4*4)* Clearly the conclusions reached 
—T- " 

regarding convergence of the sequence | are also valid for the 
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sequence | • For the sequence we are unable to make any con¬ 


clusion relative to b since the of Ca,se II is 


= a + 


TT 


—rr 

4 X 

o 


= a + 



= a + 



(4. 1.6) 


4 (a + J ) 


and lies in the interval (a,b). Moreover, we can see that for i ^ 1 
each element of the sequence is smaller than the corresponding 


element of the sequence we conclude that the limit of the se¬ 


quence 


M- 


X is smaller than x . We have that 


X = a + 


V‘+ 


= a + 


. ,a+x» 
4{^—) 


(x - a)^[ (a + x)^ + i ] - (a + x)^ = 0 



Vl + 


(a + x) 


(4. 1.7) 


2,-2 


2 ,, - 


X + (1 - 2a - ir ) X - [2a(l + it )] x + + a (l-ir ) = 0 (4. 1. 8) 

For the larger roots of (4« 1„ 8) we look at (4,, 7) and see that 

2 

(a + x) 1 so that an approximation yields 

(x - a)^ (a + x)^ - (a + x)^ ir^ = 0 
X = a + ir « 

We conclude that (4« U 8) has a solution near a + tt, and that this 
solution is the x we desireo 

Table III lists the desired roots of (4«lo8) for the same 20 roots 
as in section 4, A comparison of the error column of table III against 
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table I shows that the difference between x and the tabulated root b is 
about one-tenth as large as the difference between x and b at the first 
root. Moreover, we see that x approaches b more rapidly than x, 
that X - b is only . 01 times ^ - b for the tenth root. 

B ased upon the results of these computations, we conclude that 
although X is nearer b than x , x is still not the value of x whose 
existence is guaranteed by theorem 2. 2* 

4.2. Mean Value. 

Consider the integrated mean value of P(x) over the interval (a,x) 
and denote this integrated mean value by P(x), then 

^ X 


P(x) = 


= — f 

Ja 


(1 + 


1 


4t 


) dt 


1 


x-a 

J_ 

x-a 


Ft - — 1^ 

^ 4t ■'a 

[x-a + T ( - - r ) J 


1 


= 41 ^ ) 


= 1 + 


4ax 


We see that 


(a+x)'^ 

i V 1 

■ (a+x)^ 

a^ + 2ax + x^ ^ 4ax 
a^ - 2ax + x^ ^ 0 
(a-x)^ > 0 for all x 4 

Suppose we generate a sequence ^^x.^ 
recursion formulae: 


(4.2. 1) 


a . 

by means of the following 
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1P{X.) = i + 

' 1 4ax* 


u.** + IP(x.) u. = 0 

1 ' i' 1 


Uj(x) = k. sin V^(x.)(x-a) 


^i+1 = a + ir/ VlP(xp 


(4. 2.2) 

(4.2.3) 

(4.2.4) 
(4. 2.5) 


But equations (4*. 2„ 2) through (4, 2. 5) are the same as the set of 

2 

equations (4.4) through (4.7) except ax^ is substituted for x^ in 
(4.4). If we make this substitution throughout our preceeding considera- 
tion we see that the sequence /also converges to some x = x such 


that 


x = a + 


TT 



\/l + 


4ax 


(4ax + l)(x-a)^ - 4ax = 0 

4ax^+ (l-8a^) x^ + (4a^ - 2a - 4a tt^) x + = 0 (4. 2.9) 

Inequality (4. 2. 1) shows us that x must indeed be less than x 
but once more we are unable to state any analytical conclusions regard¬ 
ing X and b . 

Table IV contains the roots of (4. 2. 9) for the same zeros of 
J^(x) as previously computed. We see that x is smaller than b. We 
note also that jx-bj is greater than |x-bj , 

Since b has been bracketed by x and x , a linear average of the 
two is suggested. Accordingly we define 
P(xl.) ^ (P(x.) + P(xp)/2 


xL^j - a + -it/ Vp(xi7) 


(4.2. 10) 
(4.2. 11) 


Then equation (4. 2. 10) and (4. 2. 11) together with equations (4. 5) and 
(4. 6) constitute a set of recursion formulae which will generate a 
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sequence |xl. ' ^ We know the sequence xi will converge to some 

xl since it is the sum of two convergent sequences arid^ x/,- . The 
values of xl are displayed in Table We see that although the values 

of xl are nearer to b than they are not as near as x « 

o— 

Suppose we consider the average of x and x . Let us define: 

P{x2.) = (P(x.) + P(x.))/2 (4.2.12) 

^^i+1 “ a + it/ V"p(x2^ ) (4.2.13) 

Again equations {4«Z„12) and (4.2.13) together with (4.5) and (4.6) 

will generate the convergent sequence ^x2.(^ . We see from Table VI 

that x2 is nearer to b than any of the previous values generated. The 

“ -3 

maximum difference between x2 and b is less than 10 for the first 

7 

root and decreases to less than 10 for the 20th root. The difference 
decreases quite ra.pidly as is evidenced by a value less than 10"^ beyond 
the fifth root. 

4.3 Cumulative Effect. 

Let us now investigate the cumulative error for each of the sequen¬ 
ces previously discussed. The errors in this case are cumulative in 
the following sense: 

Assume that a = r^, the first root of Jq(x) is known, we will 
compute an x from a sequence. Then each successive root is generated 
using the root last computed as the value of the known root. 

The cumulative errors for x, x , x, xl and x2 are shown in 
Tables VII through Xl respectively. As well expected, the cumulative 
error is due mostly to the error in computing the first root. We see from 
Table XI tha,t the cumulative error for the 20th root is only 10% 
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greater than the error for the first root^ The worst case is shown in 
Table X where the cummulative error tor the 20th root is 40% greater 
than the error for the first root« 

* 
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(5. 1) 


5o Application to BesseT s Equation of Order n ^ 1 . 

When n J>. 1, the normal form of Bessel's equation is 

y»t + (1 + ) y = 0 

4x 

Let us impose the condition 
y{a) = y(b) = 0 . 

Again let P{x) denote the coefficient of y in (5. 1), i.e, 

1 -4n^ 

P(x) = 1 + (5. 2) 

4x __ 

We see that P(x) > 0 for all x moreover, P(x) is 

monotone increasing, lim P(x) = 1. i.e^ 0 < P(x) < 1 for all x 

X j ,— - 

such that x is in the interval V4n^-1, . We shall limit our dis¬ 

cussion to the case P(x) > 0 in order that the conditions of the theorems 
of section 2 may be satisfied^, 

Except for P(x^), the recursion formulae defining the sequence will 
be the same as in section 4, i,e. 


P(x ) = 1 + 

(5. 3) 

4x 


o 


u!' + P(x.) u. = 0 

1 ' i' 1 

(5.4) 

u^(x) = k. sin^P(x^) (x-a) 

(5.5) 

X.. , = a + ^/ Yp(x.) 

1+1 ^ ' i' 

(5.6) 


As before we consider the two cases x = a and = a + ir , 

o o 

Case I: 

Choose X = a , giving 


TD/ X _ 1 l-4n _ , ^ l~4n 

^ , z- - ^ + -~r 

4x 4a 

o 

(5. 7) 

u*' + P{x ) u =0 
o o o 

(5.8) 

u^(x) =• sinVP(x^) (x-a) 

(5.9) 

Xj = a + IT / \/P(x^) 

(5. 10) 
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Since P(x) is increasing, we see that P(x^) < P(x) for all x in 

(a,b). Then theorem 2. 1 tells us that b . Moreover, ^ > 

therefore, from (5, 10) we see that > a + tt , 

When X, is substituted in the recursion equation we have 
^ 2 

P{x ) = 1 + (5. 11) 

4Xj 

Uj* + P{xj) Uj = 0 

u^{x) = kj sin Vp( X J ) (x-a) 

= a + ir / P{xj^) (5. 12) 

We know from the inequality x^^b that x^ is in the interval (a,b). 
Trivially then < Xj, moreover, P{Xj) < 1 tells us that x^ > a + tt « 

Then we have 


X =a<a+'rr<x^<b<Xi 
o 2 1 


(5. 13) 


Let us now consider the next element of the sequence. If x^ is 


substituted in equations (5« 2) through (5, 6) ,we have 
P(x ) = 1 + 

^-2 


(5. 13) 


u 


2 + ^<^2^ ^"2 = ° 


u,(x) k, sin Vp(x,) (x-a) 


Xj = a + tr/Vpiix^) 


‘ 2 '''' ■ ■'2 ° ' 2 ' 

(5. 14) 

Since is in the interval (a,b) we do not have the necessary 

inequality for theorem 2, 1 and consequently we can no longer determine 
the relation of x^ to b for i ^ 3 , We see, however, that PCx^ )>P(Xj) 
since P(x) is increasing and < Xj . This implies that x^ ^ x^ . 

Clearly P^x^) < P{x^) and x^ < Xj „ Then we have 


X = a < a + IT < Xt < Xo < X, 
o c, .i I 


(5. 15) 
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Suppose, for seme ^^3 that the following inequalities hold: 


a < a f tr^K^ < u . o < ^ ^ j 

The recur sion formula gives 


< X 


7r~ 

k.2 

2 


P(Xk) 


1 + 


P(Xk_i) = i + 


4 K 
l,-4ri 

1 ^ ^ 
I -4n 


(5, 16) „ 

(5. 17) 

(5. 18) 

(5, 19) 


4x 


k- I 


'^k-2 ^ ^^^k-2)^k»2 " 0 

{5. 20) 

uj^ + ° 

(5.21) 

^k-1 ° 

(5.22) 

1 ” St f IT / V^P(x.^ p) (x-a,) 

(5.23) 

Xk = a + TT / Yp{Xk.^^) (x-a) 

(3. 24) 

■^k+1 ‘ ■'* ^ (x-a) 

(5. 25) 


The assumption (5„ 16) tells v* s that Pi^k 2 ^ ^PtXk)^P(Xk j) since 
P(x) is increasing. If we apply theorem 2, 1, first with respect to (5. 21) 
and (5,20), then with respect to (5„ 2i) and (5„ 22); we have 

= a < a f w < >2 < = ., < Xk_2 ^k+l^ ^k- l'^" ” “ "^^1 

We have now shown inductively that lies in the interval 


(x., X. ,) for all 1 . L.el us now select the sub sequence |z .f of lx. ,? 

ii-1 

such that for each z. - « Trivially “x „ Assuming 

j 1 o o 

for some m>0 we see that z 


then z y z 
m 


m 




^ m-1 

Mo reover^ we know that lies in the interval (x, , ,) 

V 3 mir 1--1 tj m ^ m i i 


therefore a we conclude z 


mf"! 


y ^ 


m 


and that ifz .!> is monotone increasing^ 
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Clearly, the sequexice iZj is bounded abo^/^i Dy Therefore, we 

conclude that the sequence Zj converges f.o some z „ 


Similarly, consider the subsequence of chosen such that 

T' 1 

for each n, argument similar to the one for 'we 

can show that the subsequence is decreasing and bounded below by 

a, therefore, it converges to some v « We have now shown that the 
sequence Ix/| contains two convergent subsequences JZj ^^and ;v^ ; , and 
in order to prove convergence we must show^ that the subsequences 

converge to the same point, Co that z ^ We know from (5„ 25) 
that 


a. 4’ 




a f 


TT 


1^4n 


(5. 26) 


fl + 

/ 4 / a + 


J 


V 




ir 

aT il! 

4 7 " 


and V 


a f 


/ 


I-in 


4 / , + ^ 

f 


V 




/i 1 i-4n 

' V ‘ ^ ; 


\ 


Equation (5,26). after considerable algebraic manipulation, can be 
reduced to the eqni.'Ion 

[z-a)^‘(4 |'a|4zt- i-4n' ) f 4z | + (I-4n )(4z 1-4n^')-4rr‘‘ 

|a(4z l-4n'^) f- 4? | 16a‘‘i’"(4z^ ' (z- a)“'- 4Tr'';' 0„ (p, 28) 

which is clearly a polynomial of degree 8 „ Similarly, Equation (5„ 27) . 
will reduce to the same polynomiaJ, Then we can conclude thad z - v 


and the sequence |>.^ | converge s to some, x such that 
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^ I 

X == a + 


which reduces to 


IT 


1 + 


1.4n 


4x 


4'x ^-8ax ^ + (4a^+ l«4n^->4ir^)x ^-2a(l-4n^) x 

+ (l-4n^) = 0 (5. 29) . 

Moreover, x is the root of equation {5« 29) nearest a + tt . 

Case II: 

Suppose we choose = a + ir then (5„ 3) through (5. 6) becomes 

P(x ) = 1 + = 1 + ^ (5. 28) 

^ 4x^ 4(a+7r) 

u** + P(x ) u - 0 
o o o 


( 5 ^ 29 ) 


u^(x) = sin Y P{x^) (x-a) 

X^ = a + TT / V 


We have chosen an x^ with interval (a,b) , therefore, we cannot 
make any comparisons between x^ and b. Since 
P{x ) < 1 we know x, > a + tt = x 
Let us now find the next element of the sequence . We 

have the equations 


p(-i) - ^ 


2 


4x 


1 


uy + P(vj) Uj = 0 


Uj(x) s kj sin VP(xj) (x-a) 



The inequality x, >x implies that P{x J < P(x ) since P(x) is 


s therefore, tells us that < Xj^ , 


decreasing. Theorem 2 
moreover, from P(> j| < 1 we see that x^ ^ a + tt giving 


X 


a + IT < X 2 < 


(5. 30) 
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Clearly^ ior ail i> Z this sequence hv-^Laves exactly as the sequence 


of Case I, therefore^ converges to some Xj such 


that 


Xj - a + 


ir 


1 + 


1 -4n 


4 'vy 

X 


(5.31) 

But equation (5»3^1} after simplification reduces to the same quartic as 
{5o27), and we concJude that = x « 

In either case for a > 1 we have been unable to relate x to b, 
therefore, we must consider computations for the first 20 roots* 

The results of the computation displayed in table XII are com¬ 
parable to those shown in table T for the case n s 0, as should well be 

expected* We see that for the 20th root, the error in table XII is 

-4 - 6 

10 while for J , it is about 5 x 10 
o 

Suppose, however, the sequence to converge to b* The 

final result would for some the following 

rr 

b = a f 

"P(,t } 


= a + 


1 + 


l-4n 


4 ". ' 

(AX ^ \ L 4o“) (b-a)- - A:i ^ 

which when solveo lor % snows itiat 
'U- ^ a) 



((b-a)^- ir^) 


(5. 32) 


Trivially t. dt^^irtniaed in this manner is real •■aiued, since 
(b-a)^-iT and 4n^ ^ 1 . The values of % computed from equation (5.32) 
are shown in table XIII. These values suggest that the following modifi¬ 
cation be made to the recursion formulae. 
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P(x.) 


i + 


1 4n^' 

. .a- fy, / 

4( 


u " + P{x.)a = 0 

1 S3 


u^(x) - sin YP(x^){x-a) 


X.. , s a + 
1+1 


IT 


/p(x.) 


We have already shown that the sequence ;^x/f defined by equations 

'< 3 . 


(5,3), (5,4)j (5„ 5) and (5^6) converges to x. The sequence is 

the sunn of a constant and a multiple of x^, and therefore will converge. 
This implies that the sequ^-nc e |converges to some x such that 


X = a + 


I + 


1-n 

M^) 


(5. 13) 

x^+(l-4n^-2a^-Tr^")x' - [Za|i *4n^+tT^) ]x + a'^f - ir^) - 0 (5. 14) 

The roots of the quartc (So 14) nearest a r are displayed in 
table XIV for th< first aU) roots of Jj{x)o Comparing the errors of 
table XIV against Uvble Xll we see that jx>bj is one order of magni¬ 
tude smaller than j>-bj o 

We shall 'amine the sequence | i, io^o th^ mean value of 

P(x) over the interv.^! (a^ The mean value P(y) is defined in the 
same manner as in clion 4^2^ We know that the sequence | f will 
converge since it consists of elements which are multiples of the elements 
of a convergent sequent ^ o Then we have 


= a + it/ y P(x 

ir 


a + 




1 


1^4r 

yv 

4a X 


29 







• • * 










(5. 15) 


4a^^+{l-4n^-8a^)x*"+(4a^-2a«8a^"n^-4aTT)x a - 0 

We see that for n ^ 1, the limit of the sequence results in a cubic 
similar to the one for x when n = 0, Again we know that x is the root 
of the cubic nearest a + tt. 

Table XV displays the roots of equation (5„ 15) for the first 20 
roots of Jj^{x). As expected from the results for we see that 

again x is past b and that x and x have straddled b. We note, 
as for Jq{x), Ix-bj > |x-b| . 

We shall now examine xl and x2, as defined in section 4. 2. The 

values of xl are displayed in table XVI , Again, as for J^{x), we see 

that is further from b than x« 

x7 is again nearest b for Jj(x), with a maximum error for the 

first root. In general, we see that the vaiious sequences generated for 

Jj{x) tend to converge to points farther from b than the same sequence 

for J (x). 
o 

As expected, the cumulative errors shown in tables XVIII through 
XXII exhibit the same characteristics as the cumulative errors for 
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6. Conclusions. 


We have proposed a technique of computing the consecutive zeros 
of a solution to a Sturm-Liouville System. 

[r(x)y'(x)]' + P(x)y(x) = 0 (6. 1) 

y(a) = 0 y(b) = 0 

The technique requires the determination of an ^ such that the equation 

[r(x)u'(x)]' + P(x)y(x) = 0 (6.2) 

will have a solution u{x) which maybe forced to vanish at x = a and 

at X = b . 

The existence of such an x has been shown by theorem 2. 2. More¬ 
over, in theorem 2. 2« 1 we have established upper and lower bounds on 
the range of this x . 

In the specific case of the Bessel functions* we note that each of the 

sequences generated by the proposed techniqu-^ co iverges to a particular 

X, no matter whether the first element was chosen equal to the known 

root a, or chosen equal to a + ir . We see that x turns out to be the 

root of a quartic nearest a + tt. Similarly^we see tha.t x is also a root 

/\ 

of a quartic polynomial while x resolves into a cubic equation. 

The solution for % for the case n = 0 and for n = 1 both indicate 
that the value of rhf: desired x as defined by theorem 2^ 2 is almost 
midway between a and bo The difference between x and X approaches 
Tr/2 as x incr }ases. 

We have also shov/n that the desired x for the Bessel functions lies 
between the x determined by a linear averaging technique and the x 
determined by the mean value of P(x) over the interval (a,x)o 

Unfortunately, all of the above generating techniques failed to result 
in a sequence converging to the desired x. 
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TABLE 1 


VALUES OF X *** Jq(X) 


ZEROS 

OF THE BESSEL FUNCTION OF ORDER ZERO 

FOR ROOTS 

ONE THROUGH TWENTY, BASED ON THE RECURSION 

FORMULAE 

P(X(n) = 1 + 1/4(X(I > 

X(I+1> » 

A+TT/ \/p(X( 1) ) 

ERROR 

IS THE DIFFERENCE 

BETWEEN TABULATED 

AND COMPUTED ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

5.5200781103 

5.5336719233 

.0135938131 

2 

8.6537279126 

8.6564432415 

.0027153292 

3 

1.1 .7915344390 

11.7925004759 

.0009660370 

4 

14.9309177082 

14.9313671612 

.0004494532 

5 

18.0710639674 

18.0713085639 

.0002445965 

6 

21.2116366294 

21.2117842031 

.0001475737 

7 

24.3524715304 

24.3525673207 

.0000957905 

8 

27.4934791317 

27.4935447974' 

.0000656658 

9 

30.6346064676 

30.6346534281 

.0000469606 

10 

33.7758202134 

33.7758549489 

.0000347356 

n 

36.9170983527 

36.9171247659 

.0000264135 

12 

40.0584257636 

40.0584463133 

.0000205500 

13 

43.1997917127 

43.1998080136 

.0000163012 

14 

46.3411883712 

46.3412015187 ' 

.0000131483 

15 

49.4826098960 

49.4826206546 

.0000107595 

16 

52.6240518400 

52.6240607537 

.0000089143 

17 

55.7655107528 

55.7655182229 

.0000074702 

18 

58.9069839241 

58.9069902431 

.0000063197 

19 

62.0484691886 

62.0484745828 

.0000053946 

20 

65.1899647992 

65.1899694391 

.0000046415 
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table II 

VALUES OF **• Jq(X) 


ZEROS 

OF THE BESSEL FUNCTION OF ORDER ZERO 

FOR ROOTS 

ONE THROUGH TWENTY, BASED 

ON THE RECURSION 

FORMULAE 

CHI 

IS THE VALUE X MUST 

ASSUME TO CONVERGE TO 8 

5G = ((8- 

•A)/2)( Vi/((B-A)^< 

-7i2) ) 


ERROR 

IS THE DIFFERENCE 

BETWEEN X AND CHI 


ROOT 

NO, 


TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 


5.5200781103 

3.8368678503 

1.6968040731 

2 


8.6537279126 

7.0180586502 

1.6383845913 

3 


11.7915344390 

10.1751194785 

1.6173809974 

4 


14.9309177082 

13. 3249362609 

1.6064309003 

5 


18.0710639674 

16.4716300704 

1.5996784936 

6 


21.2116366294 

19,6166964089 

1,3950877941 

7 


24.3524715304 

22.7608070695 

1.5917602512 

8 


27.4934791311 

25.9043037482 

1.5892410493 

9 


30.6346064676 

29.0473744823 

1.5872789458 

10 


33.7758202134 

32.1902150838 

1.5856398651 

11 


36.9170983527 

35.3327245601 

1.5844002059 

12 


40.0584257636 

38.4751889762 

1.5832573371 

. 13 


43.1997917127 

41.6175294435 

1.5822785702 

14 


46.3411833712 

44.7596824104 

1.5815191083 

15 


49,4826098960 

47.9016739726 

1.5809466820 

16 


52.6240518400 

51.0439078268 

1.5801529270 

17 


55.7655107528 

54.1856175233 

1.5799006997 

18 


58.9069839241 

57.3278564233 

1.5791338198 

19 


62.0484691886 

60.4698329102 

1,5786416726 

20 


65.1899647992 

63.6115127150 

1.5784567241 
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TABLE III 

VALUES OF Z *** Jq(X) 


ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMU LAE 

p(x(i)) = i+i/u( ( (A+jTd ) )/2)? x(i + i) * A+Tt/ y/?iZii)) 


ERROR 

IS THE DIFFERENCE 

BETWEEN TABULATED 

AND COMPUTED ROOTS 

ROOT 

NO. 

TABULATED 

ROCT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5217120022 

.0016338920 

2 

8.6537279126 

8.6538810786 

.0001531661 

3 

11.7915344390 

11.7915695005 

.0000350617 

4 

14.9309177082 

14.9309296827 

.0000119747 

5 

18.0710639674 

18.0710691069 

.0000051396 

6 

21.2116366294 

21.2116391868 

.0000025579 

7 

24.3524715304 

24.352472942B 

.0000014123 

8 

27.4934791317 

27.4934799736 

.0000008423 

9 

30.6346064676 

30.6346070007 

.0000005335 

10 

33.7758202134 

33.7758205654 

..0000003527 

11 

36.9170983527 

36.9170985967 

.0000002443 

12 

40.0584257636 

40.0584259368 

.0000001732 

13 

43.1997917127 

43.1997918384 

.0000001260 

. 14 

46.3411883712 

46.3411884652 

.0000000946 

15 

49.4826098960 

49.4826099686 

.0000000730 

16 

52.6240516400 

52.6240518950 

.0000000556 

17 

55.7655107528 

55.7655107975 

.0000000452 

18 

58.9069839241 

58.9069839586 

.0000000349 

19 

62.0484691886 

62.0484692166 

.0000000280 

20 

65.1899647992 

. 65.1899648216 

.0000000226 
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table IV 

VALUES OF $ *** Jq{X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTYt BASED ON THE RECURSION FORMULA 
P($(I)J = UI/4A)^(n ^(I + l) » A+T/ 

ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 


I- 


x( n) 


ROOT 

NO* 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5172325156 

-.0028455946 

2 

8.6537279126 

8.6534819095 

-.0002460030 

3 

11.7915344390 

11.7914791510 

-.0000552879 

4 

14.9309177082 

14.9308989581 

-.0000187501 

5 

18.0710639674 

18.0710559469 

-.0000080201 

6 

21.2116366294 

21.2116326448 

-.0000039842 

7 

24.3524715304 

24.3524693325 

-.0000021976 

8 

27.49347913lt 

27.4934778223 

-.0000013094 

9 

30.6346064676 

30.6346056401 

-.0000008275 

10 

33.7758202134 

33.7758196630 

-.0000005496 

n 

36.9170983527 

36.9170979746 

-.0000003774 

12 

40.0584257636 

40.0584254945 

-.0000002689 

/ 13 

43.1997917127 

43.1997915152 

-.0000001969 

14 

46.3411883712 

46.3411882240 

-.0000001468 

15 

49.4326098960 

49.4826097042 

-.0000001110 

16 

52.6240518400 

52.6240517525 

-.0000000871 

17 

55.7655107528 

55.7655106848 

-.0000000672 

18 

58.9069839241. 

58.9069838692 

-.0.000000548 

19 

62.0484691886 

62.0484691439 

-.0000000444 

20 

65.1899647992 

65.1899647620 

-.0000000366 
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TABLE V 

VALUES OF in *•* Jq{X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE_ 

pixTd)) = (P(x(n)+p(${i)))/2 irr(i+i)» A+^ 'vpTxTiiT) 

ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 


ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5254224626 

.0053443524 

2 

8.6537279126 

8.6549613329 

^0012334206 

3 

11.7915344390 

11.7919896499 

.0004552110 

u 

14.9309177082 

14.9311330232 

.0002153153 

5 

18.0710639674 

18.0711822445 

.0001182773 

6 

21.2116366294 

21.2117084200 

.0000717907 

7 

24.3524715304 

24.3525183247 

.0000467948 

8 

27.4934791317 

27.4935113089 

.0000321774 

9 

30.6346064676 

30.6346295336 

.0000230662 

10 

33.7758202134 

33.7758373059 

.0000170928 

11 

36.9170983527 

36.9171113698 

.0000130179 

12 

40.0584257636 

40.0584359039 

.0000101405 

13 

43.1997917127 

43. 1997997640 

•0000080521 

14 

46.3411883712 

46.3411948718 

.0000065007 

15 

49.4826098960 

49,4826152194 

.0000053242 

16 

52.6240518400 

52.6240562536 

.0000044135 

17 

55.7651^107528 

55.7655144539 

.0000037015 

.18 

58.9069839241 

58.9069870561 

.0000031325 

19 

62.0484691886 

62.0484718634 

.0000026751 

20 

65.1899647992 

65.1899671014 

.0000023024 
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TABLE VI 


VALUES OF X2 *** JqIX) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 

ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE_ 

pilcTd)) =(p(7( I) )+p{x( I)) )/2 «(i + i) » A+'^/VpufTnT” 

ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 


ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5194687883 

”.0006093219 

2 

8.6537279126 

8.6536814596 

-.0000464527 

3 

11.7915344390 

11.7915243239 

-.0000101150 

U 

14.9309177082 

14.9309143201 

-.0000033879 

5 

18.0710639674 

18. 0710625267 

-.0000014403 

6 

21.2116366294 

21.2116359160 

-.0000007132 

7 

24.3524715304 

24.3524711374 

-.0000003927 

8 

27.4934791317 

27.4934788980 

-.0000002336 

9 

30.63U6064676 

30.6346063204 

-.0000001470 

10 

33.7758202134 

33.7758201146 

-.0000000984 

1 1 

36.9170983527 

36.9170982856 

-.0000000665 

12 

40.0584257636 

40.0584257152 

-.0000000478 

13 

43.1997917127 

43.1997916764 

-.0000000355 

14 

46.3411883712 

46.3411883442 

-.0000000261 

15 

49.4826098960 

49.4826098764 

-.0000000190 

16 

52.6240518400 

52.6240518242 

-.0000000158 

17 

55.7655107528 

55.7655107416 

-.0000000110 

18 

58.9069839241 

58.91369839139 

-.0000000099 

19 

62.0484691886 

62.0484691802 

-.0000000082 

20 

65.1899647992 

65.1899647918 

-.0000000070 
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TABLE VII 

CUMULATIVE ERRORS FOR “x *** Jq(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER the first ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 
FROM THE VALUE CF 'x LAST OBTAINED. ERROR IS CUMULATIVE 
FROM THE FIRST ROOT. 


ROOT 

NO* 

TABULATED 

ROCT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5336719233 

.0135938131 

2 

8.6537279126 

8.6700534131 

*0163255008 

3 

11.7915344390 

11.8088337616 

.0172993229 

4 

14.9309177082 

14.9486705526 

.0177528444 

5 

18.0710639674 

18.0390637650 

«0179997976 

6 

21.2116366294 

21.229785478? 

.0181488496 

7 

24.3524715304 

24. 3707171553 

.0182456251 

8 

27.4934791317 

27.5117911105 

.0183119791 

9 

30.6346064676 

30.6529659065 

.0183594390 

10 

33.7758202134 

33.7942147609 

.0183945483 

11 

36.9170983527 

36.9355196003 

.0184212479 

12 

40.0584257636 

40.0768677853 

.0184420225 

13 

43.1997917127 

43.2182502151 

.0184585024 

14 

46*341iCe3712 

46.3596601672 

.0184717962 

15 

49.4826098960 

49.5010925708 

.0184826752 

16 

52.6240518400 

52.6425435282 

.0184916886 

17 

55.7655107528 

55.7840099940 

.0184992421 

18 

58.9069839241 

58.9254895560 

.0185056320 

19 

62.0484691886 

62.0669802753 

,0185110874 

20 

65.1899647992 

65.2084805798 

.0185157806 
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TABLE VIII 

CUMULATIVE ERRORS FOR X *** Jq(X) 


ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 


AFTER 

THE FIRST ROOT 

IS COMPUTED, THE NEXT 

ROOT IS COMPUTED 

FROM 

THE VALUE OF X 

LAST OBTAINED. ERROR 

IS CUMULATIVE 

FROM THE FIRST ROOT. 



ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

5.520078 1103 

5.5217120022 

.0016338920 

2 

8.6537279126 

8.6555185516 

.0017906391 

3 

11.791534L390 

11.7933614515 

.0018270126 

l4 

14.9309177082 

14.9327572957 

.0018395875 

5 

18.0710639674 

18. 0729090152 

.0018450481 

6 

21.2116366294 

21.2134844260 

.0018477967 

7 

24.3524715304 

24.3543208619 

.0018493315 

8 

27.4934791317 

27.4953293884 

.0018502570 

9 

30.6346064676 

30.6364573166 

.0018508493 

10 

33.7758202134 

33.7776714582 

.0018512453 

11 

36.9170983527 

36.9189498741 

.0018515220 

12 

40.0584257636 

40.0602774834 

.0018517202 

13 

43.1997917127 

43.2016435778 

.0018518660 

14 

46.3411883712 

46.3430403462 

.0018519759 

15 

49.4826098960 

49.4844619567 

.0018520612 

16 

52.6240518400 

52.6259039668 

.0018521272 

17 

55.7655107528 

55.7673629336 

.0018521811 

18 

58.9069839241 

58.9088361468 

.0018522234 

19 

62.0484691886 

62.0503214458 

.0018522573 

20 

65.1899647992 

65.1918170843 , 

.0018522854 
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TABLE IX 

CUMULATIVE ERRORS FOR X *•* Jq(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 
FROM THE VALUE OF $ LAST OBTAINED. ERROR IS CUMULATIVE 
FROM THE FIRST ROOT. 


ROOT 

NO. 

TABULATED 

ROCT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5172325156 

-.0028455946 

2 

8.6537279126 

8.6506294168 

-.0030984957 

3 

11.7915344390 

1 1.7883782729 

-.0031561660 

U 

14.9309177082 

14.9277417252 

-.0031759829 

5 

18.0710639674 

18.067879399S 

-.0031845679 

6 

21.2116366294 

21.2084477427 

-.0031888865 

7 

24.3524715304 

24.3492602321 

-.0031912980 

8 

27.4934791317 

27.4902863791 

-.0031927526 

9 

30.6346064676 

30.6314127846 

-.0031936828 

10 

33.7758202134 

33.7726259045 

-.0031943080 

n 

36.9170983527 

36.9139036089 

-.0031947432 

12 

40.0584257636 

40.0552307060 

-.0031950568 

13 

43.1997917127 

43. 1965964232 

-.0031952894 

14 

46.3411883712 

46.3379929066 

-.0031954644 

15 

49.4826098960 

49,4794142973 

-.0031955985 

- 

52.6240516400 

52.6208561352 

-.0031957048 

17 

55.7655107528 

55.7623149641 

-.0031957879 

16 

58.9069339241 

58.9037880665 

-.0031958566 

19 

62.0484691886 

62.0452732751 

-.0031959134 

20 

65.1899647992 

65.1867688391 

-.0031959599 
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table X 

CUMULATIVE ERRORS FOR xT *** Jq(X) 


ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 


AFTER THE FIRST ROOT 

IS COMPUTED. THE NEXT 

ROOT IS COMPUTED 

FROM 

THE VALUE OF xi 

LAST obtained. ERROR 

IS CUMULATIVE 

FROM 

THE FIRST ROOT. 



ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

5.5200781103 

5.5254224626 

.0053443524 

.2 

8.6537279126 

8.6603153765 

.0065874640 

3 

11.7915344390 

11.7985812158 

.0070467769 

4 

14.9309177082 

14.9381818194 

.0072641112 

5 

18.0710639674 

18.0784474835 

.0073835164 

6 • 

21.2116366294 

21. 2190926266 

.0074559972 

7 

24.3524715304 

24.3599747745 

.0075032441 

8 

27.4934791317 

27.5010148650 

.0075357334 

9 

30.6346064676 

30.6421654909 

.0075590234 

10 

33.7758202134 

33.7833964955 

.0075762823 

11 

36.9170983527 

36.9246877786 

.0075894268 

12 

40.0584257636 

40. 0660254285- 

.0075996651 

13 

43.1997917127 

43.2073995080 

.0076077954 

14 

46.3411883712 

46.3488027304 

.0076143594 

15 

49.4826098960 

49.4902296308 

.0076197354 

16 

52.6240518400 

52.6316760313 

.0076241914 

17 

55.7655107528 

55. 7731386814 

.0076279288 

18 

58.9069839241 

58.9146150155 

.0076310916 

19 

62.0484691886 

62.0561029809 

.0076337926 

20 

65.1899647992 

65. 1976009160 

.0076361171 
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TABLE XI 

CUMULATIVE ERRORS FOR 1(2 *** Jq ( X ) 


ZEROS 

OF THE BESSEL FUNCTION OF ORDER ZERO 

FpR ROOTS 

ONE THROUGH TWENTY, BASED 

ON THE RECURSION 

FORMULAE 

AFTER THE FIRST ROOT IS 

COMPUTED, THE NEXT ROOT IS COMPUTED 

FROM 

THE VALUE OF X2 LAST OBTAINED. ERROR 

IS CUMULATIVE 

FROM 

THE FIRST ROOT. 



ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

5.5200781103 

5.5194687883 

-.0006093219 

2 

8.6537279126 

8.6530707316 

-.0006571810 

3 

11.7915344390 

n. 7908666495 

-.0006677894 

4 

14.9309177082 

14.9302463082 

-.0006714000 

5 

18.0710639674 

18.0703910086 

-.0006729586 

6 

21.2116366294 

21.2109628869 

-.0006737421 

7 

24.3524715304 

24.3517973502 

-.0006741801 

8 

27.4934791317 

27.4928046870 

-.0006744443 

9 

30.6346064676 

30.6339318538 

-.0006746134 

10 

33.7758202134 

33.7751454851 

-.0006747281 

11 

36.9170983527 

36.9164235452 

-.0006748069 

12 

40.0584257636 

40.0577508984 

-.0006748646 

13 

43.1997917127 

43.1991168046 

-.0006749080 

14 

46.3411883712 

46.3405134305 

-.0006749401 

15 

49.4826098960 

49.4819349311 

-.0006749644 

16 

■ 52.6240518400 

52.6233768547 

-.0006749846 

17 

55.7655107528 

55.7648357525 

-.0006749997 

10 

58.9069839241 

58.9063089108 

-.0006750130 

19 

62.0484691886 

62.0477941642 

-.0006750239 

20 

65.1899647992 

65.1892897654 

-.0006750335 
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TABLE XII 

VALUES OF *X *** (X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE through twenty, BASED ON THE RECURSION FORMULAE_ 

p(x(i)) = i+(i-4N^)/u(x(i ))^ x(i + i) = a+71/ /pTxTm 


ERROR 

IS THE DIFFERENCE 

BETWEEN TABULATED 

AND COMPUTED ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

7.0155866698 

6.9976376109 

•“.0179490588 

2 

10.1734681348 

10.1686351451 

-.0048329895 

3 

13.3236919360 

13.3217202593 

-.0019716767 

4 

16.4706300506 

16.4696368496 

-.0009932009 

5 

19.6158585101 

19.6152890925 

-.0005694173 

6 

22.7600853802 

22.7597279320 

-.0003574477 

7 

25.9036720875 

25.9034342724 

-.0002378150 

8 

29.0468285335 

29.0466620042 

-.0001665289 

9 

32.1896799095 

32.1895587798 

-.0001211293 

10 

35.3323075492 

35.3322166996 

T.0000908489 

n 

38.4747662339 

38.4746963540 

-.0000698795 

12 

41.6170942122 

41.6170393117 

-.0000549003 

13 

44.7593189972 

44.7592750816 

-.0000439151 

14 

47.9014608869 

47.9014252098 

-.0000356767 

15 

51.0435351832 

51.0435058065 

-.0000293767 

16 

54.1855536401 

54. 1855291631 

-.0000244770 

17 

57.3275254359 

57.3275048267 

-.0000206090 

16 

60.4694578443 

60.4694403261 

-.0000175173 

19 

63.6113566970 

63.6113416851 

-.0000150119 

20 

66.7532267310 

66.7532137688 

-.0000129620 
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table XIII 

VALUES OF 9^ *** Jj{X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 


CHI IS THE VA LUE X MUST ASSUME TO CONVERGE TO B 
% - ( {B-A)/2) { >/{ l-4N^)/( (B-A ) 


ERROR 

IS THE DIFFERENCE 

BETWEEN X AND CHI 


ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROCT 

ROOT 


1 

7.0155866698 

5.3312873839 

1.6663502270 

2 

10.1734681348 

8.5374782400 

1.6311569051 

3 

13.3236919360 

1 1.7071 112888 

1.6146089705 

4 

16.4706300506 

14. 8645464401 

1.6050904095 

5 

19.6158585101 

18.0163550223 

1.5989340702 

6 

22.7600853802 

21. 1610798943 

1.5986480376 

7 

25.9036720875 

24.3119656080 

1.5914686644 

8 

29.0468285335 

27.4576333771 

1.5890286271 

9 

32.1896799095 

30.6024408266 

1.5871179532 

10 

35.3323075492 

33. 7466444736 

1.5855722260 

n 

38.4747662339 

36.8904200587 

1.5842762953 

12 

41.6170942122 

40.0338375270 

1.5832017846 

13 

44.7593189972 

43. 1769980257 

1.5822770558 

14 

47.9014608869 

46.3199358247 

1.5814893851 

15 

51.0435351832 

49.4627229534 

1.5807828531 

16 

54.1855536401 

52. 6053818800 

1.5801472832 

17 

57.3275254359 

55.7479370888 

1.5795677379 

18 

60.4694578443 

58. 890202 9432 

1.5792373829 

19 

63.6113566970 

62.0326487795 

1.5786929056 

20 

66.7532267310 

65. 1751772780 

1.5780364908 
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TABLE XIV 
VALUES OF r *** 


Jj^(X) 


ZEROS OF THE BESSEL FUNCTION 
ONE THROUGH TWENTY, BASED ON 
PIXdJ) = 1 + (1-UN^)/4((A+X(I) 
ERROR IS THE DIFFERENCE BETWE 


OF ORDER ONE FOR ROOTS NUMBER 
THE RECURSION FORMULAE 
)/2)^ 7ll + n = A+7Z/ /pixi I ) ) 

EN TABULATED AND COMPUTED ROOTS 


ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

7.0155866698 

7.0141418219 

-.0014448478 

2 

10.1734681348 

10.1732513588 

-.0002167758 

3 

13.3236919360 

13.3236309085 

-.0000610273 

U 

16.4706300506 

16. 4706066190 

-.0000234314 

5 

19.6158585101 

19.6158476649 

-.0000108451 

6 

22.7600853802 

22.7600786900 

-.0000066901 

7 

25.9036720875 

25.9036688204 

-.0000032669 

8 

29.0468285335 

29.0468265270 

-.0000020061 

9 

32.1896799095 

32.1896786084 

-.0000013006 

10 

35.3323075492 

35.3323066691 

-.0000008800 

1 1 

38.4747662339 

38.4747656174 

-.0000006158 

12 

41.6170942122 

41.6170937680 

-.0000004442 

13 

44.7593189972 

44.7593186684 

-.0000003282 

14 

'47.9014608869 

47.9014606383 

-.0000002479 

15 

51.0435351832 

51.0435349923 

-.0000001903 

16 

54.1855536401 

54.1855534911 

-.0000001483 

17 

57.3275254359 

57. 3275253186 

-.0000001170 

18 

60.4694578443 

60.4694577483 

-.0000000957 

19 

63.6113566970 

63.6113566197 

-.0000000766 

20 

66.7532267310 

66.7532266676 

-.0000000615 
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table XV 

VALUES OF X *** J_^(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 

ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE^_ 

P{5(I)) » 1 ♦ { 1 - 4N^)/ UAX( n ^(I + l) » A I )) 

ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 


ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

7.0155866698 

7.0180468041 

.0024601344 

2 

10.1734681348 

10.1738162260 

.0003480912 

3 

13.3236919360 

13.3237883153 

.0000963795 

u 

16.4706300506 

16.4706667955 

.0000367450 

5 

19.6158535101 

19.6158754556 

.0000169458 

6 

22.7600853802 

22.7600932522 

.0000078722 

7 

25.9036720875 

25.9036771739 

.0000050868 

8 

29.0468285335 

29.0468316572 

.0000031241 

9 

32.1896799095 

32.1896819314 

.0000020222 

10 

35.3323075492 

35.3323089145 

.0000013660 

n 

38.4747662339 

38.4747671895 

.0000009564 

12 

41.6170942122 

41.6170949005 

.0000006892 

13 

44.7593189972 

44.7593195057 

.0000005093 

14 

47.9014608869 

47.9014612706 

.0000003842 

15 

51.0433351832 

51.0435354784 

.0000002958 

16 

54.1855536401 

54.1855538711 

.0000002316 

17 

57,3375254359 

57.3275256194 

.0000001642 

18 

60.4694576443 

60.4694579896 

.0000001462 

19 

• 63.6113566970 

63.6113568163 

.0000001198 

20 

66.7532267310 

•66. 7532.268297 

,0000000997 






TABLE XVI 

VALUES OF ri *** J|(X) 


ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FqR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULA^_ 

pcriti)) - (p( x( I) )+p (^( I)))/? jTfd) « M+7C/V piTTi 1 )) 


ERROR 

IS THE DIFFERENCE 

BETWEEN TABULATED 

AND COMPUTED ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

7.0155866698 

7.0077901660 

-.0077965038 

2 

10.17346815411 

10.1712216930 

-.0022464416 

3 

13.3236919360 

13.3227535994 

-.0009383365 

4 

16.470630C506 

16.4701516442 

-.0004784064 

5 

19.6158585101 

19.6155822147 

-.0002762953 

,6 

22.7600853802 

22.7599105686 

-.0001748114 

7 

25.9036720875 

25.9035557127 

-.0001163747 

8 

29.0468285335 

29.0467468258 

-.0000817076 

9 

32.1896799095 

32.1896203524 

-.0000595564 

10 

35.3323075492 

35.3322628057 

-.0000447430 

11 

38.4747662339 

38.4747317713 

-.0000344625 

12 

41.6170942122 

41.6170671051 

-.0000271062 

13 

44.7593189972 

44.7592972936 

-.0000217033 

14 

47.9014608869 

47.9014432402 

-.0000176465 

15 

51.0435351832 

51.0435206424 

-.0000145406 

16 

54.1855536401 

54.1855415171 

-.0000121228 

17 

57.3275254359 

57. 3275152231 

-.0000102125 

18 

60.4694578443 

60.4694491578 

-.0000086856 

19 

63.6113566970 

63.6113492502 

-.0000074461 

20 

66.7532267310 

66.7532202993 

-.0000064312 
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TABLE XVU 

VALUES OF ^2 ♦** J2^{X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE_ 

= (p(x( I >)+p(^(i)) )/2 )^(i) = a+TT/ \/p c iTz 11 )) 


ERROR 

IS THE DIFFERENCE 

BETWEEN TABULATED 

AND COMPUTED ROOTS 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

7.0155866698 

7.0160913781 

.0005047083 

2 

10.1734681348 

10.1735337211 

.0000655864 

3 

13.3236919360 

13.3237096060 

,0000176702 

4 

16.4706300506 

16.4706367063 

.0000066559 

5 

19.6158585101 

19.6158615602 

.0000030501 

6 

22.7600853802 

22.7600859711 

.0000005910 

7 

25.9036720875 

25.9036729969 

,0000009099 

8 

29.0468285335 

29.0468290923 

.0000005590 

9 

32.1896799095 

32.1896802699 

.0000003608 

10 

35.3323075492 

35.3323077913 

.0000002430 

11 

38.4747662339 

38.4747664034 

.0000001703 

12 

41.6170942122 

41.6170943342 

.0000001225 

13 

44.7593189972 

44.7593190875 

.0000000905 

14 

47.9014608869 

47.9014609549 

.0000000681 

15 

51.0435351832 

51.0435352353 

.0000000527 

16 

54.1855536401 

54.1855536811 

.0000000416 

17 

57.3275254359 

57.3275254695 

.0000000336 

18 

60.4694578443 

60.4694578694 

.0000000253 

19 

63.6113566970 

63.6113567185 

.0000000216 

20 

66.7532267310 

66.7532267496 

.0000000191 
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TABLE XVIII 

CUMULATIVE ERRORS FOR 'J? »** JjIX) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 


FROM 

THE VALUE OF 

LAST OBTAINED. ERROR 

IS CUMULATIVE 

FROM 

THE FIRST RCOT. 



ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

7.0155866698 

6.9976376109 

“.0179490588 

2 

10.1734681348 

10.1507267661 

-.0227413687 

3 

13.3236919360 

13.2990017352 

-.0246902007 

4 

16.4706300506 

16.4449597476 

-.0256703027 

5 

19.6158585101 

19.5896268408 

-.0262316693 

6 

22.7600853802 

22.7335015247 

-.0265838551 

7 

25.9036720875 

25.8768550311 

-.0268170560 

8 

29.0468285335 

29.0198475332 

-.0269810002 

9 

32.1896799095 

32.1625796901 

-.0271002193 

10 

35.3323075492 

35.3051179312 

-.0271896177 

11 

38.4747662339 

38.4475078629 

-.0272583706 

12 

41.6170942122 

41.5897818329 

-.0273123789 

. 13 

44.7593189972 

44.7319634212 

-.0273555757 

14 

47.9014608869 

47.8740702206 

-.0273906654 

15 

51.0455351832 

51. 0161 156254 

-.0274195571 

16 

54.1855536401 

54.1581100114 

-.0274436282 

17 

57.32?5254359 

57. 3000615416 

-.0274638942 

18 

60.4694578443 

60.4419767251 

-.0274811186 

19 

63.6113566970 

63.5838608174 

-.0274958793 

20 

66.7532267310 

66.7257181071 

-.0275086237 
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table XIX 

CUMULATIVE ERRORS FOR X **• JiIX) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 


FROM 

THE VALUE OF X 

LAST OBTAINED. ERROR 

IS CUMULATIVE 

FROM 

THE FIRST RCOT. 



ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 

1 

7.0155866698 

7.0141418219 

-.0014448478 

2 

10.1734681348 

10. 1718119476 

-.0016561872 

3 

13.3236919360 

13. 3219771462 

-.0017147896 

4 

16.4706300506 

16.4688930572 

-.0017369931 

5 

19.6158585101 

19.6141 113704 

-.0017471395 

6 

22.7600853802 

22.7583319838 

-.0017533960 

7 

25.9036720875 

25.9019167111 

-.0017553763 

8 

29.0468285335 

29.0450713504 

-.0017571828 

9 

32.1896799095 

32.1879215697 

-.0017583393 

10 

35.3323075492 

35.3305484364 

-.0017591120 

11 

38.4747662339 

38. 4730065875 

-.0017596461 

12 

41.6170942122 

41.6153341858- 

-.0017600260 

13 

44.7593189972 

44.7575586932 

-.0017603031 

14 

47.9014608869 

47.8997003762 

-.0017605102 

15 

51.0435351832 

51.0417745160 

-.0017606668 

16 

54.1855536401 

54.1837928528 

-.0017607870 

17 

57.32/5254359 

57.3257645555 

-.0017608804 

IS 

60.4694578443 

60.4676968884 

-.0017609558 

19 

63.6113566970 

63.6095956815 

-.0017610151 

20 

66.7532267310 

66.7514656689 

-.0017610620 





TABLE XX 

CUMULATIVE ERRORS FOR ^ *** JjIX) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE through twenty, BASED ON THE RECURSION FORMULAE 


AFTER 

THE FIRST ROOT 

IS COMPUTED, THE NEXT 

ROOT IS COMPUTED 

FROM 

THE VALUE OF ^ 

LAST OBTAINED. ERROR 

IS CUMULATIVE 

FROM THE FIRST ROOT. 



ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

7.0155866698 

7.0180468041 

.0024601344 

2 

10.1734681348 

10.1762664455 

.0027983109 

3 

13.3236919360 

13.3265823792 

.0028904433 

U 

16.4706300506 

16.4735551225 

.0029250719 

5 

19.6158585101 

19.6187993325 

.0029408227 

6 

22.7600853802 

22.7630333365 

.0029479564 

7 

25.9036720875 

25.9066256429 

.0029535557 

8 

29.0468285335 

29.0497848745 

.0029563413 

9 

32.1896799095 

32. 1926380284 

.0029581190 

10 

35.3323075492 

35.3352668518 

.0029593028 

1 1 

38.4747662339 

38.4777263533 

.0029601197 

12 

41.6170942122 

41.6200549109 

.0029606996 

13 

44.7593189972 

44.7622801177 

.0029611211 

14 

■'47.9014608869 

47.9044223214 

.0029614345 

15 

51.0435351832 

51.0464968551 

.0029616724 

16 

54.1855536401 

54.1885154955 

.0029618555 

17 

57.3275254359 

57.3304874348 

.0029619992 

18 

60.4694578443 

60.4724199548 

.0029621108 

19 

63.6113566970 

63.6143188979 

.0029622011 

20 

66.7532267310 

66.7561890054 

.0029622754 


I 
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TABLE XXI 

CUMULATIVE ERRORS FOR 7 \ *** J]^(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FQR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 


FROM 

FROM 

THE VALUE OF XI 

THE FIRST RCOT. 

LAST OBTAINED. ERROR 

IS CUMULATIVE 

ROOT 

TABULATED 

COMPUTED 

ERROR 

NO. 

ROOT 

ROOT 


1 

7.0155866698 

7.0077901660 

“.0077965038 

2 

10.1734681348 

10.1634497307 

“.0100184040 

3 

13.3236919360 

13.3127478317 

-.0109441041 

U 

16.4706300506 

16.4592144503 

-.0114155999 

5 

19.6158585101 

19.6041707369 

-.0116877730 

6 

22.7600853802 

22.7482254347 

“.0118599452 

7 

25.9036720875 

25.8916985528 

-.0119735343 

8 

29.0468285335 

29.0347745535 

-.0120539797 

9 

32.1896799095 

32. 1775672976 

-.0121126117 

10 

35.3323075492 

35. 3201508904 

-.0121566580 

11 

38.4747662339 

38.4625756508 

-.0121905831 

12 

41.6170942122 

41.6048769467 

-.0122172653 

13 

44.7593189972 

44.7470803680 

-.0122386289 

lU 

47.9014608869 

47. 8892048877 

-.0122559989 

15 

51.0435351832 

51.0312648714 

-.0122703116 

16 

54.1855536401 

54.1732713953 

-.0122822442 

17 

57.3275254359 

57.3152331384 

-.0122922968 

18 

60.4694578443 

60.4571569972 

-.0123008466 

19 

63.6113566970 

63.5990485204 

-.0123081761 

20 

66.7532267310 

66.7409122232 

-.0123145066 
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TABLE XXII 

CUMULATIVE ERRORS FOR )(2 **• JjIX) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 
FROM THE VALUE OF Zz LAST OBTAINED. ERROR IS CUMULATIVE 
FROM THE FIRST RCOT. 


ROOT 

NO. 

TABULATED 

ROOT 

COMPUTED 

ROOT 

ERROR 

1 

7.0155866698 

7.0160913781 

.0005047083 

2 

10.1734681348 

10.1740364628 

.0005683281 

3 

13.3236919360 

13.3242770866 

.0005851508 

4 

16.4706300506 

16.4712214330 

.0005913828 

5 

19.6158585101 

19.6164527028 

.0005941927 

6 

22.7600853802 

22.7606000152 

.0005946353 

7 

25.9036720875 

25.9042685339 

.0005964467 

8 

29.0468285335 

29.0474254708 

.0005969373 

9 

32.1596799095 

32.1902771583 

.0005972489 

10 

35.3323075492 

35.3329050038 

.0005974551 

11 

38.4747662339 

38.4753638301 

.0005975968 

12 

41.6170942122 

41.6176919080 

.0005976967 

13 

44.7593189972 

44.7599167656 

.0005977688 

14 

47.9 u.*.605 86 9 

47.9020587085 

.0005978224 

15 

51.0435351832 

51.0441330457 

.0005978626 

16 

54.1855536401 

54.1861515343 

.0005978945 

|7 

57.0275254359 

57.3281233553 

.0005979196 

18 

60.4694578443 

60.4700557813 

.0005979377 

19 

63.6115566970 

63.61 19546490 

.0005979527 

20 

66.7532267310 

66.7538246959 

.0005979659 
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